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We show that the difference resonance driven by the space charge pseudo-octupole of high- intensity 
beams not only couples the beam core emittances; it can also lead to emittance exchange in the 
beam halo, which is of relevance for beam loss in high intensity accelerators. With reference to 
linear accelerators the "main resonance" k z /k XtV = 1 (corresponding to the Montague resonance 
2Q X — 2Q y — in circular accelerators) may lead to such a coupling and transfer of halo between 
planes. Coupling of transverse halo into the longitudinal plane - or vice versa - can occur even if the 
core (rms) emittances are exactly or nearly equal. This halo argument justifies additional caution 
in linac design including consideration of avoiding an equipartitioned design. At the same time, 
however, this mechanism may also qualify as active dynamical halo cleaning scheme by coupling a 
halo from the longitudinal plane into the transverse plane, where local scraping is accessible. We 
present semi-analytical emittance coupling rates and show that previously developed linac stability 
charts for the core can be extended - using the longitudinal to transverse halo emittance ratio - to 
indicate additional regions where halo coupling could be of importance. 



I. INTRODUCTION 

The possibility of manipulating the 6D phase space 
distribution of particle beams by exchanging emittances 
between planes can be important for optimum perfor- 
mance as well as for minimizing beam loss. In the context 
of free electron lasers, for example, emittance exchange 
between transverse and longitudinal planes has been pro- 
posed to take advantage of very small source longitudi- 
nal emittances. For this purpose dispersion-based con- 
cepts for multi-GcV electron beams involving chicanes, 
transverse RF cavities as well as quadrupolcs have been 
suggested [l[ . For (partially stripped) ion beams the pro- 
posal was made recently that a beam with equal trans- 
verse emittances can be effectively transformed into a 
more flat beam matched better to subsequent ring injec- 
tion, where this rotator includes normal quadrupoles, a 
stripper, chicane and skew quadrupoles [2J. 

For circular machines the idea of employing lattice non- 
linear resonances to work on beam tails and halo clean- 
ing was proposed some time ago With reference to 
the CERN Intersecting Storage Ring the suggestion was 
made to use a 5 th order resonance to "pump" protons 
into the tail up to the aperture or a distant scraper. This 
would allow to clean the halo dynamically rather than by 
applying - as usual - a scraper close to the beam core. 
In a different context coupling between tails was em- 
ployed experimentally at the CERN Proton Synchrotron, 
where insufficient (transverse) Landau damping of head- 
tail modes in one plane could be successfully enhanced 
by linear coupling with the tail distribution in the other 
plane plane 

In high intensity accelerators any active concept ad- 
dressing emittance exchange needs to deal with space 
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charge and possible side effects of it. In our study we con- 
sider the intrinsic space charge pseudo-octupole - always 
present in a non-uniform spatial distribution - as driv- 
ing term of a nonlinear difference resonance and explore 
its possible role in an active emittance coupling scheme. 
Core and halo emittance coupling is an important issue 
for high intensity linacs, where trade-offs between longi- 
tudinal acceptance and synchronous phase require careful 
consideration (see, for example, Ref. [5[). 

In linear accelerators this space charge difference res- 
onance is known as "main resonance" with the equiva- 
lent resonance condition in linac notation k z /k x ^ y = 1. 
For a discussion of its theory see Refs. [H-IH; furthermore 
Ref. 9] for the first experimental observation of the linac 
main resonance at the GSI UNILAC. For reference we 
mention that in circular accelerators this mechanism is 
better known as "Montague resonance" [l(| with the cor- 
responding difference resonance condition 2Q X —2Q y = 0. 
For a detailed experimental study of the Montague reso- 
nance in the CERN Proton Synchrotron see Ref. [n| . 

In discussions of high intensity linac beams it is often 
assumed that resonant interaction with a space charge 
driven nonlinearity primarily affects the beam core. For 
the latter it is assumed that rms emittance transfer 
between longitudinal and transverse planes should be 
avoided. This can be achieved either by using a pri- 
ori "equipartitioned" beams, or simply by avoiding the 
corresponding resonance condition and work with non- 
equipartitioned beams and a higher flexibility of lattice 
tuning. The latter option has found new interest for H~- 
linacs, where avoiding the risk of intrabeam stripping 
losses [HI is helped by reduced transverse focusing and 
removing the constraint of an equipartition requirement. 

The main point of this study is to explore whether 
the emittance coupling concepts describing primarily the 
beam core can be extended to describe coupling effects in 
the beam halo - whatever the origin of it may be. Along 
this line we also examine the possibility of employing 
space charge and the core-halo coupling in a strategy for 
active transfer of longitudinal halo into the transverse 
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plane, where it can be scraped in a well-controlled way. 

Our study is carried out with the TRACEWIN code 
in an idealized transport lattice with RF focusing, hence 
ignoring acceleration and all additional complexity of a 
specific linac lattice. In Section [TT] we present some de- 
tails of the ideal model; in Section HiTl results for dynam- 
ical crossing of the main resonance with/ without halo; 
in Section IIVI we present results for fixed tunes; in Sec- 
tion [V] analytical scalings are summarized; in Section WT\ 
we introduce extended stability charts, with concluding 
remarks in Section [VIII 



II. LATTICE, BEAM CORE AND HALO 
DEFINITIONS 

Without lack of generality we employ an "idealized ac- 
celerator" lattice with a periodic FODO channel (85 cells, 
with cell length arbitrarily chosen as 1 m) with linearly 
varying transverse phase advance and - for simplicity - a 
single RF kick per cell to generate longitudinal focusing 
with a fixed phase advance of 74° per cell throughout the 
paper (Fig. []] only showing first six cells) . The current 
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FIG. 1. TRACEWIN envelopes for reference transport chan- 
nel with RF focusing gaps. 

is chosen such as to yield a transverse tune depression 
of about 20%, which is modest for high current linacs. 
We note that for this kind of space charge physics ac- 
celeration can be ignored, because the only parameters 
that matter are the ratio of phase advances per cell, the 
tune depression and the ratios of normalized core or halo 
omittances. Since we always assume that the initial omit- 
tances and phase advances in x and y are equal, we only 
indicate ratios of z and x quantities at the start. 

We make the following model assumptions: (i) The ini- 
tial distribution follows the TRACEWIN option of ran- 
domly distributed particles in a six-dimensional ellipse 
(here truncated at 2.8 a); (ii) 10 5 simulation particles and 
a 3D Poisson solver (PICNIC); (hi) an "extra halo" dis- 
tribution is introduced as six-dimensional Gaussian type 
cut at 2 a and - if present - occupied by 2.5% of the total 



number of particles; (iv) the halo size is expressed by the 
99.9% emittance as area of an ellipse defined by using 
core twiss parameters and containing 99.9% of the parti- 
cles; (v) we allow for a nonsymmetric (anisotropic) halo 
with independent emittances transversely and longitudi- 
nally. The halo emittances are described by independent 
multiplication factors M z /M xy with respect to the core 
emittances. The value of 2.5% halo intensity is relatively 
arbitrary and higher than what can usually be tolerated 
in high power linacs. Yet it is still low enough to ensure 
that the halo space charge force is only a weak contribu- 
tion, which is the main argument for it in our study. 

Note that this study assumes initial Gaussian halo pro- 
files throughout and ignores entirely the question of what 
physical mechanism generates the halo, which might for 
example be mismatch or errors in linear accelerators. A 
projection of a characteristic initial core plus halo distri- 
bution into the real plane z — x is shown in Fig. [5] for 
M z /M x = 6/1.5. 
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FIG. 2. Exemplary initial core-halo distribution in z — x plane 
for halo multiplication factors M z /M x = 1.5/6, with ellipse 
including 99.9% of total particles. 



III. SIMULATION RESULTS FOR DYNAMICAL 
CROSSING 

In this section we present simulation results under the 
assumption of a linearly decreasing transverse phase ad- 
vance such that the k z /k x = 1 is crossed from left to 
right. We mention that the analogous crossing in the 
opposite direction yields similar but not equal results. 
There are some subtle differences, however, which will 
be discussed further in Section [TVl 
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A. Crossing without extra halo 

For reference we first review this crossing in the ab- 
sence of extra halo as shown in Fig. [5] (M z /M xy = 1/1). 
The tune footprint with initial rms emittances e z /e x = 
1.5 is plotted on the stability chart in the top of Fig. [3] In 
this chart for an initial rms emittance ratio e z /e x = 1.5 
the stop-band increases with space charge, which pro- 
vides the driving term. For a detailed discussion see Sec- 
tion EJ. The footprint is determined in TRACEWIN as 
ratio of the k z and k xy shown in the insert in the upper 
left corner of Fig. [3] these values are extracted as rms 
values of tunes obtained from tracking of all particles. A 
linear tune sweep k xy o : 85° — > 70° over 85 cells (here 
85 m) and initial tune depression k x jk§ x = 0.8 are as- 
sumed. The crossing of the center of resonance condition 
k z /k xy = 1 in this case occurs at cell 50. Note that in 
this and all following graphs the distance on the abscissa 
in meters is identical with the cell number for our lattice 
example. 

It is noted that the crossing leads to the expected 
rms emittance transfer, where the final state is nearly 
"equipartitioned" . The longitudinal rms emittance re- 
duction is about double the increase of each transverse 
one. This reflects the fact that the longitudinal degree of 
freedom has to "heat" both transverse - a kind of energy 
conservation supported by the observation that the sum 
of all three rms emittances is constant within few per- 
cent. The initial values of the 99.9% emittances reflect 
a small halo due to truncation of the core distribution 
at 2.8 a. They follow the same pattern of exchange as 
the rms emittances, including an approximate constancy 
of their sum. For the following discussion it is impor- 
tant to note that the equipartition process itself does not 
generate additional halo - only exchange. 
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FIG. 3. Crossing of main resonance from left to right in the 
absence of extra halo: Tune footprint on stability chart for 
£z/€x = 1.5 (top, with graphs for k z and k xy in the upper left 
corner)), Rms (center) and 99.9% (bottom) emittances. 



B. Extra purely transverse halo 

Next we consider the same parameters for the core as 
before, but now we add a purely transverse halo distribu- 
tion with multiplication factors M z /M xy — 1/4.5, which 
results in an initial ratio e z /e x — 0.4 for the 99.9% emit- 
tances and inverts the core and halo emittance ratios. 
The tune ramp is extended to k X Q = 85° — > 65° over 
again the 85 cells. The rms emittances again roughly 
equipartition. This also applies to the 99.9% emittances, 
but their transfer occurs into the opposite direction: the 
two initially larger transverse halo emittances "pump" 
the longitudinal halo emittance, which more than doubles 
in a first swing. Thereafter, apparently, a small fraction 
of halo particles oscillates back and forth between planes 
- approximately around halo emittance equipartition as 
shown in Fig. 2J 



C. Nearly equipartitioned core and purely 
longitudinal halo 

Here we assume an almost equipartitioned core with 
tz/^x — 0.95, which results in a practically negligible 
stopband width for the main resonance as shown in 
Fig. [5j We furthermore add a purely longitudinal halo 
with corresponding factors M z /M xy = 4.5/1, which im- 
plies a starting ratio e z /e x = 3 for the 99.9% emittances. 
The transverse tune is ramped from k X Q — 85° — >• 70 . 
The coupling has a weak effect on the rms emittances, 
which merge together as well - partly due to the shrink- 
ing longitudinal halo. The 99.5% emittances come to 
a full exchange with some overshoot and qualitatively 
the same behavior as the rms and 99.5% emittances in 
Fig. [3] Hence, the longitudinal halo is almost fully con- 
verted into a transverse halo with little effect on the core 
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(top), 99.9% (bottom) emittances (center) and projection of '■ 
z-evolution (bottom). 20 40 60 80 
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emittances. Fig. [5] also suggests that the coupling of halo 
emittances is subject to a much broader stopband than 
the one shown for the core, which is calculated for the 
core emittance ratio (see also Section IVII 



FIG. 5. Extra halo in longitudinal plane with nearly equipar- 
titioned core: Fune footprint on stability chart for e z /e x = 
0.95 (top); Rms and 99.9% emittances (middle) and z- 
evolution (bottom). 



D. Dynamical halo cleaning with transverse 
scrapers 

Here we assume an isotropic halo with M z /M xy = 
4.5/4.5. As expected no emittance exchange occurs, as 
this beam is equipartitioned both in core and halo. We 
next define a transverse scraper by introducing in three 
consecutive cells circular apertures at positions, where 
the beam is round. In our example such a scraper of 
4 mm radius sharply eliminates the transverse halo par- 
ticles (with 3% reduction of beam intensity). Assum- 
ing that the scraper is located in the lattice before the 
coupling resonance is reached, the latter can pump the 
longitudinal halo into the transverse plane. The loca- 
tion of this scraper is found to be most effective at cells 
38/39/40 (equal to distance in m), where k z /k x = 0.93. 
Note that the tune ratio reaches the center of the stop- 



band (k z /k x = 1) after slightly less than 20 additional 
cells. For our parameters 20 cells is also the characteris- 
tic emittance exchange time for a static tune placed right 
on the stopband, hence without ramping across it. The 
issue of characteristic emittance exchange time is further 
discussed in Section IIVI as it is the basis for choosing the 
scraper position most effectively. 

The result of this simulation is shown in Fig. [51 It is 
seen that the coupling of the longitudinal halo into the 
transverse plane makes the transverse 99.9% emittances 
grow - following again the already observed conservation 
of the sum of all three emittances. To avoid a partial re- 
turn into the longitudinal plane, which would bring the 
99.9% e z back to about the value 1.5 in our example, we 
have added a second scraper with only two apertures in 
cells 75/76. The second scraper makes sure that the halo 
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cleaning remains conserved for all three planes. The rms 
emittances are only changing by a few percent as a result 
of the scraping. Note that the second scraper is unneces- 
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FIG. 6. Cleaning of initially isotropic halo by two consecu- 
tive transverse scrapers: 99.9% emittances (top); x-evolution 
(center) and z-evolution (bottom). 

sary, if fine tuning of the tune ramp is realized such that 
the coupling stopband is left after the longitudinal halo 
emittance has reached its minimum. This helps reducing 
the scraped intensity - also in view of the fact that some 
remaining transverse halo is not an issue of big concern 
in linacs with large enough transverse aperture. Remov- 
ing the longitudinal halo this way appears promising as 
it helps shifting the synchronous phase towards higher 
acceleration rate. 



IV. SIMULATION OF STATIC TUNES 

In this section we study the topology of core and halo 
coupling at the main resonance stopband for fixed values 
of k xy o ("static tunes"), which reveals some interesting 
details not resolvable with dynamical crossing. In par- 
ticular, the validity of different stopband widths for core 



and halo is clearly seen in cases where their emittance 
ratios differ. 



A. Same emittance ratio in core and halo 

Keeping the transverse tunes fixed reveals different fea- 
tures of exchange, depending on the distance from the 
stopband and at which side. We assume e z /e x = 3 and 
M z /M xy — 4.5/4.5, which gives the halo the same ini- 
tial emittance ratio of 3:1. With k xy o = 75° the cor- 
responding k z /k xy starts at the right edge of the stop- 
band as shown in Fig. [7] The starting rms emittance 
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FIG. 7. Selfconsistent evolution with k xy o = 75 : tune foot- 
print on (zoomed) stability chart for e z /e x = 3 (top, with 
graphs for k z and k XtV in the upper right corner); Rms (cen- 
ter) and 99.9% (bottom) emittances. 

exchange leads to a shrinking e z /e x , which leads to corre- 
sponding changes in space charge densities with the effect 
that k z /k xy also shrinks and moves further left into the 
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stopband; this self-consistent mechanism continues until 
equality of emittances occurs - with even some overshoot 
- and k z /k xy reaches the left edge of the stopband and 
returns back to the center (as indicated by an arrow in 
the top graph of Fig. [7]). The 99.9% emittances, however, 
show only weak coupling. 

A very different behavior is found for k xy o = 70°, for 
which k z /k xy is actually outside the right edge of the 
stopband as shown in Fig. [8] Not surprisingly, there is 
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FIG. 8. Selfconsistent evolution with k xy o = 70°: tune foot- 
print on stability chart for e z /e x = 3 (top, with graphs for fc z 
and k XtV in the upper left corner); Rms (center) and 99.9% 
(bottom) emittances. 

practically no rms emittance exchange; however, the halo 
emittances show effective although not complete cou- 
pling. To interpret this phenomenon we note that tunes 
of halo particles differ from core particle tunes: the tune 
ratio of very large amplitude halo particles approaches 
asymptotically the space-charge- free value k z o/k xy Q. On 
the other hand, the corresponding resonance condition 



main resonance stopband, which helps understand the 
dominant halo particle coupling there. In fact, we have 
not found any effect on halo coupling for a working point 
at or beyond the l.h.s. edge of the main resonance stop- 
band in Fig. HI 

The asymmetric response on the main resonance stop- 
band is also found in the following: for k z /k xy starting 
close to the left edge of the stopband the core emittance 
exchange shifts k z /k xy outside of the stopband and the 
process stops - in contrast with the behavior at the right 
edge as in Fig. [71 Note that for an initial longitudinal 
emittance smaller than the transverse one the side be- 
havior is also reversed. 



B. Symmetric core with longitudinal halo 

Assuming e z /e x = 0.95 we expect practically no ef- 
fect on the rms emittances at or near k z /k xy — 1 as 
already found in the crossing case of Fig. [5j A purely 
longitudinal halo with M z /M xy = 4.5/1 shows emittance 
exchange depending on the value of k z /k xy . We find that 
halo omittance coupling is most effective at the condition 
k z o/k x o = 1 (see Fig. [5]) and gradually vanishes away 
from it in either direction. Note that this is supported 
by the fact that for equal rms emittances the condition 
k z /k xy = 1 coincides with the condition k z o/k X Q = 1 - 
the asymptotic resonance condition for very large ampli- 
tude halo particles - in contrast with the case of unequal 
emittances in Section HVAI 
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FIG. 9. Static tune with k z o/k x o 
(bottom) emittances. 



1: Rms (top) and 99.9% 
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By testing different vaiues of k z /k xy we actually find 
that the region of halo emittance exchange is described 
quite well by a stability chart generated with the halo 
emittance ratio (rather than core emittance ratio). This 
can be explained by the observation that the driving term 
for both the exchange of core and halo particles is the 
space charge multipole. This finding will be pursued fur- 
ther in Sections |V] and |VI] 



V. SEMI-ANALYTICAL SCALINGS 

For practical applications it may be useful to consider 
the approximate scaling expressions for the main reso- 
nance stopband width and exchange time, which have 
been derived in Ref. Using Ak = ko — k and the 
tune ratio k z /k x as variable, the stopband width can be 
written in the form 



(1) 



where it is noted that the longitudinal tune depression is 
(approximately) related to the transverse one according 
to 




Ak z /k z0 



Ak x /k xa 



(2) 



As an example we mention that for e z /e x — 1.5 and 
k x /k x o — 0.6 the resulting S fts 0.18 is in good agreement 
with Fig. H 

A corresponding expression for the e-folding time for 
emittance coupling is suggested in Ref. [7J for the center 
of the stop-band. By comparison with simulation results 
we find that in the region of moderate emittance imbal- 
ance 0.5 < e z /e x < 2 the number of cells required to reach 
the first crossing of emittances is basically independent 
of the emittance ratio and found as 



N ceUs « - (Ak z /k z0 ) 1 . 



(3) 



We also find that the coupling distance estimates in Eq. [3] 
apply quite well to the halo exchange in spite of the fact 
that the coupling of halo particles with the core - via its 
space charge pseudo-octupole - is expected to decrease 
with increasing distance from the core. 



VI. EXTENDED STABILITY CHARTS 



individual particles without a net global transfer. It is 
thus appropriate to use distinct stability charts for core 
and halo, where the rms emittance ratio is used for the 
core, and the 99.9% (or similar) emittance ratio for the 
halo. This is demonstrated in Fig.[T0]for a crossing of the 
main resonance from left to right with k xy o : 85° — > 70°; 
furthermore the assumption of a nearly equipartitioned 
core (e z /e x — 0.95) and a purely transverse halo with 
M z /M xy = 1/4.5, which results in a 99.9% emittance ra- 
tio e z /e x = 3.7. We define an extended stability chart in 
the following way: the halo stability chart is plotted as 
contour lines on top of the solid core chart - each one with 
the corresponding emittance ratio. The 99.9% emittance 
curve gives clear evidence that this exchange occurs over 
a broad range of k z /k xy , which we find is comparable 
with the width of the stopband (halo contours) at the 
main resonance. 
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FIG. 10. Crossing of main resonance: tune footprint on ex- 
tended stability chart combined for core (e z /e x = 0.95) and 
halo (e z /e x — 3.7)(top) and 99.9% emittances (bottom). 



Our calculations give clear evidence that halo emit- 
tance exchange may occur independent of whether a core 
omittance exchange takes place or not. The conditions 
for effective halo emittance exchange at the main reso- 
nance are: (i) the presence of a space charge octupole as 
driving term (always given for a non-uniform beam); (ii) 
fulfilment of the resonance condition and (iii) different 
emittances, otherwise there is exchange of amplitudes of 



Crossing from right to left with k xy o : 65° — > 85° con- 
firms a similar width of response in the variable k z /k xy , 
which is shown in Fig. [TTJ Hence in both cases the very 
narrow main resonance stopband of the core is found ir- 
relevant for the halo. For creating these extended charts 
we have assumed that the above described behavior at 
the main resonance is also applicable in an analogous 
way to all other resonances on the chart. 
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FIG. 11. Same as Fig. [10] with opposite crossing of main 
resonance. 



VII. CONCLUSION 

We have shown that the presence of the space charge 
pseudo-octupole - in all non-uniform beam models - is not 
only a source of emittance exchange for the core emit- 
tances, but also for halo emittances. The strength of 
this driving term is considerable - only very few betatron 
periods are needed in a high current linac to reach full 
exchange of emittances. 

These findings allow to extend the interpretation of 
linac stability charts in order to describe separately pos- 
sible core and halo emittance exchange. For the main 
resonance the corresponding width is readily estimated 
by an approximate expression containing the respective 
emittance ratio. 

For practical linac design there is an incentive to avoid 
undesirable transfer of halo from the transverse into the 
longitudinal, which can get lost out of the bucket. This 
enhances the importance of avoiding this coupling reso- 
nance - or even the vicinity of it - in the main part of 
a linac. We find no argument at this level in support of 
the frequently made assumption that it is advantageous 
to make the core equipartitioned. On the contrary, it 
may be appropriate to relax conditions and enable a non- 
equipartitoned bunch in the interest of strictly avoiding 
the resonance and halo coupling. Refraining from en- 
forced equipartitioning may in consequence turn out ad- 
vantageous for the design. 

On the other hand, it is suggested to consider this res- 
onance for an active scheme of longitudinal halo transfer 
into the transverse direction, where it can be scraped. In 
practice such use of the resonance should be limited to 
a transition between linac structures and avoided else- 
where. In any case our study makes it clear that besides 
halo intensity it is also important to study halo emit- 
tances independently for the longitudinal and transverse 
planes. 
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